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We investigate the topological properties of the QCD vacuum with 4 flavours of dynamical staggered fermions 
at finite temperature. To calculate the topological susceptibility we use the field-theoretical method. As in the 
quenched case, a sharp drop is observed for the topological susceptibility across the phase transition. 


1. INTRODUCTION 

The topological susceptibility is defined as 

X = J d^xd^d,{T{K^{x)K,{0))) ( 1 ) 

where K^{x) is 
9 ^ 

Kfiyx) = t57]-2 ^ 

A“ [dpAl - (2) 

and Q{x) = daK^{x) is the density of topological 
charge. Fq. (|^) uniquely defines the prescription 
for the singularity of the time ordered product 
at X —> 0 0. Determining y around in full 
QCD is an important piece of information to test 
models of QCD vacuum 0. 

1.1. The Simulation 

We have simulated the theory on a 16'^ x 4 lat¬ 
tice with four flavours of staggered fermions with 
bare mass am = 0.05. We have used the HMC 
algorithm for the updating. To be sure that topol¬ 
ogy is well thermalized we have checked the topo¬ 
logical charge of our sample of configurations by 
cooling. They are well decorrelated. In Fig. 
we display the distribution of topological charge 
of our configurations at f3 = Q/g'^ = 5.04. This 
histogram confirms that at this mass value the 
topology is well sampled 0. We have learned 
that topology is easily decorrelated when we work 
with large quark masses [0|^] . In our case also the 
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lattice spacing is large, so that updating the in- 
stanton content is faster. 

With our bare fermion mass and lattice size 
the deconfining transition is known to occur at 
Pc = 5.04 0]. We have checked this number by 
computing the Polyakov loop and the chiral con¬ 
densate as shown in Fig. 0. 



Figure 1. Distribution of topological charge Ql 
at P = 5.04. 


The topological susceptibility was measured at 
/3=5.00, 5.02, 5.04, 5.05, 5.06, 5.10. 

1.2. The Operators 

We have measured the topological susceptibil¬ 
ity by using the 2-smeared lattice topological 
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Table 1 

Topological susceptibility and TjT^ vs. (3. 


/? 

T/Tc 

x/MeV^ 

5.00 

0.9677 

1.61(43) 

5.02 

0.9804 

1.13(30) 

5.04 

1.0000 

1.21(30) 

5.05 

1.0126 

2.89(1.46) 

5.06 

1.0274 

1.63(1.22) 

5.10 

1.1110 

0.43(1.19) 
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Figure 2. Polyakov loop (squares, right scale) and 
chiral condensate (circles, left scale) across the 
phase transition on a 16^ x 4 lattice. 


charge Ql ^ 

= ~ 00.^9 (n^^^(x)npj(a:)) (3) 

and the field-theoretical method |Q| to extract x 
from the lattice susceptibility 

XL = ^ = Z^a\ + M (4) 

where Ql is the lattice total charge Ql = 
^ space-time volume. The 
multiplicative and additive renormalizations Z 
and M have been evaluated by use of the heating 
method (|,^. 

2. DETERMINATION OF a, Z AND M 

The lattice spacing was extracted by measur¬ 
ing the string tension at 13 = 5.00, 5.04 and 5.10 


on a 16'^ lattice. For the other values of (3 it 
was determined by a splines interpolation. Wil¬ 
son loops were evaluated by using smeared spatial 
links. From the value at /3c {a = 0.300(20)fm) we 
infer the critical temperature to be Tc = 164(11) 
MeV. 

The heating method to evaluate Z and M 
yields a non-perturbative determination of these 
renormalization constants HJ^. To calculate Z a 
local updating algorithm is applied on a config¬ 
uration containing a charge -1-1 classical instan- 
ton. These updatings, being local, thermalize the 
short distance fluctuations, responsible for the 
renormalization effects, and due to the slowing 
down leave large structures, like instantons, un¬ 
changed. The measurement oi Ql on such up¬ 
dated configurations yields Z Q. As Q is known, 
one can extract Z. Notice that this procedure is 
equivalent to imposing the continuum value for 
the 1-instanton charge (in the MS scheme it is 
4-1) and extracting the finite multiplicative renor¬ 
malization Z by evaluating a matrix element of 
Ql- 

The additive renormalization M is obtained in 
a similar way. We apply a few heating steps with 
a local updating algorithm on a zero-field con¬ 
figuration. Then the topological susceptibility is 
calculated. This provides the value of M, if no in¬ 
stantons have been created during the few updat¬ 
ing steps. This method agrees with the treatment 
of the singularity of Eq. (|^) . 

As explained in [|| , cooling tests must be done 
to check that the background topological charge 
has not been changed during the local heating 
(it is -1-1 in the calculation of Z and 0 in the 
evaluation of M). 
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3. RESULTS 


The values for the topological susceptibility 
and the temperature are given in table ^ In Fig. 
we show the normalized topological susceptibil¬ 
ity x(T)/x(r = 0) as a function of the tempera¬ 
ture. The value at zero temperature x(T = 0) has 
been obtained by averaging the values at T < Tc- 
In the same figure the results for the quenched 
case are shown for comparison. The sig¬ 

nal for X drops strongly when crossing the tran¬ 
sition temperature. Analogous results have been 
reported for two flavours and by using the cooling 
method 0. 

From our data we cannot yet conclude that 
the drop in presence of fermions is steeper than 
for the quenched case, especially after considering 
that our simulations have been done at low beta 
values, where a poor scaling is expected, with a 
possible resulting large systematic effect on the 
scale T/Tc. 

Comparing the full QCD results to the quen¬ 
ched case is an important issue for testing instan- 
ton liquid models of the vacuum. In view of that 
we are adding determinations on a 32^ x 8 lattice 
where we expect better control on scaling. 



Figure 3. Behaviour of the topological suscepti¬ 
bility as a function of the normalized temperature 

T/n. 
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